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FIG. 1. Instantaneous relative vorticity field of the two- layer 
model on a spherical geodesic grid with 163,842 cells. 

Climate varies on decadal and longer time scales, yet 
climate simulations advance at intervals measured in 
minutes or even seconds. In 1967, however, Edward 
Lorenz observed: "More than any other theoretical pro- 
cedure, numerical integration is also subject to the crit- 
icism that it yields little insight into the problem. The 
computed numbers are not only processed like data but 
they look like data, and a study of them may be no more 
enlightening than a study of real meteorological observa- 
tions. An alternative procedure which does not suffer this 
disadvantage consists of deriving a new system of equa- 
tions whose unknowns are the statistics themselves." [1] 
We show here that such direct statistical simulation 
(DSS) is able to reproduce statistics obtained by the tra- 
ditional route of time averaging direct numerical simula- 
tions (DNS). 

DSS of damped and driven systems may be formu- 
lated by a Reynolds decomposition of the dynamical vari- 
ables qi into the sum of a mean value and a fluctua- 
tion (or eddy): qi = (qi) + q[ with (q[) = 0. The 
first two equal-time cumulants are then q = (qi) and 
Cij = (q[q'j) = (q%qj) — (Qi)(Qj)- The equation of motion 
for the cumulants form an infinite hierarchy that must be 
truncated. The simplest closure, CE2, is to set the third 
and higher cumulants equal to zero, amounting to the 
neglect of eddy-eddy interactions. Here we extend the 
work of Ref. [2 to a two-layer quasigeostrophic model 
originally formulated on the /3-plane[3] but generalized 
to the rotating sphere. DNS is carried out in real space 
on a spherical geodesic grid. [4] A typical instantaneous 
flow is shown in Fig. [T] 

Fig. [2] compares the two-point vorticity correlation 
function as calculated by DNS and CE2. Storm tracks are 
immediately apparent from the strong mid-latitude cor- 
relations and the jets fluctuate mainly at zonal wavenum- 
ber 5 and 6. CE2 captures long-range correlations ("tele- 
connection patterns") and inter- hemisphere correlations 
are weak, as expected. A quantitative comparison of the 
zonal velocities, as calculated by DNS and CE2, is pre- 
sented in Fig. [3] where the midlatitude westerlies and the 
tropical trade winds are evident. The comparison sug- 
gests that the large-scale flows are not so nonlinear as 
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FIG. 2. The two-point second cumulant of the relative vor- 
ticity field with respect to a reference field point at latitude 
45° and along the prime meridian. Top: DNS. Bottom: CE2. 
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FIG. 3. Mean zonal velocity, averaged over both layers, as 
calculated by DNS and CE2 at two different resolutions. 

to preclude their description by the cumulant expansion. 
Stratification and shearing act together to weaken the 
nonlinearities in large-scale flows. [5] It is these features 
of the general circulation that distinguishes it from the 
highly nonlinear problem of 3D isotropic and homoge- 
nous turbulence, and makes DSS possible. 

I thank P. Kushner, T. Schneider and S. Tobias. This 
work was supported in part by NSF grant Nos. DMR- 
0605619 and PHY-0551164. 


[1] E. N. Lorenz, The Nature and Theory of the General Cir- 
culation of the Atmosphere. (WMO, Geneva, 1967). 

[2] J. B. Marston, E. Conover, and T. Schneider, J. Atmos. 
Sci. 65, 1955 (2008). 

[3] S. Lee and S. Feldstein, J. Atmos. Sci. 53, 2232 (1996). 

[4] R. Heikes and D. A. Randall, Mon. Weather Rev. 123, 
1881 (1995). 

[5] P. A. O 'Gorman and T. Schneider, Geophys. Res. Lett. 
34, 524 (2007). 


